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; Abstract 
^ , 

■ Results of the computer simulation of dilepton production from the 
, expanding pion gas created in ultrarelativistic Pb+Pb 160 GeV/n col- 

I I lisions are presented. Finite size of the expanding pion gas influences 

^ ' invariant mass spectrum of dileptons. Sensitivity of the shape of dilep- 

J;^ ■ ton mass spectrum to the initial stage of the pion gas is discussed. 

I Second order asymmetry in azimuthal distribution of dileptons is pre- 

Q^ ' dieted for non-central collisions. Non-equilibrium nature of the pion gas 

■ expansion is found to exhibit itself in rapidity distribution of produced 
, dilepton pairs. 

■ 

: 1 Introduction 

^ : 

C ■ Goal of heavy ion collision (HIC) experiments is to reveal properties of the 

>• . compressed hot nuclear matter created during the collision of heavy nuclei. 

^ I Final momentum distributions of most abundant particles - hadrons are how- 

^ \ ever influenced (if not determined) during the dilute and late freeze-out stage 

■ - - ' of the heavy ion collision. Interesting information about the early stage of 

the collision is in this case hidden by subsequent collective effects of strongly 
interacting hadrons. 

Fortunately other particles - leptons and photons produced in all stages of 
HIC process, are not influenced by the later hadronic phase. Thus distributions 
of leptons and photons can provide us with more or less direct information 
about the early stages of the collision process. 
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Recent experimental results on invariant mass distribution of dilepton pairs 
|jl|, H have attracted increased attention due to the observed excess of dilepton 
yield at regions 1.5 < M^^ < 2.5GeV (NA50) and 0.2 < M^^ < 1.5 GeV 
(CERES). Theoretical analysis on dilepton pair production P] leads to the 
conclusion that unconventional production mechanisms (e.g. modified p-meson 
mass 0, 3 or in- medium pion dispersion relation [^,0]) are necessary to explain 
present experimental data. 

In this contribution specific results of the computer simulation of dilepton 
pairs production from expanding pion gas are presented. Besides the invari- 
ant mass distribution of the produced dileptons also other properties of the 
dileptonic signal produced from the pion gas via vr+vr" annihilation channel 
are analyzed. 

Short description of the computer simulation is given in next section. In- 
fluence of the finite size of pion gas (and its initial state) on invariant mass 
spectrum of dileptons is discussed in Section 3 using results of the computer 
simulation. In Section 4 rapidity distribution of dilepton pairs obtained from 
the simulation is presented. Non-equilibrium features of the pion gas expan- 
sion are discussed in this context. Second-order azimuthal asymmetry in the 
transversal momentum distribution of dileptons is predicted by the simula- 
tion of non-central Pb-|-Pb 160 GeV/n collisions in Section 5. Possible conse- 
quences of the effects found and presented here on other areas of HIC research 
are scetched at the end of this paper. 



2 Simulation of Dilepton Production 

For the simulation of dilepton production from the expanding pion gas rescat- 
tering model ^ (inspired by the model of T.J.Humanic [|10]) was used. 
Initial stage of pion gas in Pb+Pb 160 GeV/n collisions (positions, momenta 
and time of creation of pions) was taken from events generated by cascade 



program [|lT|. In the rescattering program pions move in small time steps 
AT = O.lfm as free relativistic particles and their momenta are changed in 
mutual collisions. In each time step distance of pions dij = \xi — Xj\ is com- 



pared to critical distance dc{s) = ^a.„{s)/'K, where cr^(s) is total elastic tttt 
cross section. If condition dij < dc{s) is fulfilled colhsion procedure for a given 
pair of pions is called. More detailed description of the simulation can be found 
in 1,0. 

Two independent methods of the implementation of dilepton production 
via 7r"'"'7r~ annihilation channel were used in the simulation. First one is based 
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on standard expression for the rate of dilepton production |T3[]: 



where o"e(M) is the cross section for tt+tt^ — > e+e^ at given ^/s = M, s = 



{ki + k2y, Vr is the relative velocity Vr = \j s{s — Am^) /2E1E2 and f{k,x) is 
distribution function of pions in phase space, x = (x, t). Distribution function 
f{k, x) is calculated by rescattering program in our case. It depends mainly on 
the isospin averaged cross section for elastic vrvr interaction which regulates 
scattering of pions and it is influenced also by formation time Tf and interaction 
time ti parameters of the rescattering model used 0. 

After multiplying and dividing right-hand side of Eq.(|l]) by ^^^(M^) and 
integrating over d^x the rate of dileptons from the pion gas can be expressed 
as H: 

dNee ^ <Te{M^) dN^ 

dM^ ' (T^(M2) ' dM^ ^ ^ 

where is distribution of tttt collisions as a function of invariant mass M of 
the colliding tttt pairs and R denotes ratio of the number of 7r"^7r~ collisions to 
all TTTT interactions. Distribution is calculated by the rescattering program. 
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Fig.l Distribution of pions in the reaction plane before and after the rescat- 
tering process in CMS frame of the collision. Axis values are in fm units. 

At the same time another method for the simulation of tt+tt^ annihilation in 
the expanding pion gas was used. Cross section o"e(s) for the process tt+tt^ — 
e~^e~ was multiplied by factor 13000 in the rescattering program and pairs of 
pions were tested (in each time step) for the critical distance dl{s) calculated 
from this enlarged cross section: 



dl{s) = a/13000 ■ cxe{s)/n > |x, - Xj\ (3) 
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When this condition is fulfilled for a given pair of pions invariant mass M 
and momenta of pions PiiPj are remembered (stored into file). Positions and 
momenta of pions fulfilling condition (j^) are left unchanged in the program. In 
this way probability for the process of pion annihilation is enhanced to the level 
comparable with the elastic pion-pion interactions. At the same time process 
of the expansion of pion gas is not infiuenced by this high probability of pion- 
pion annihilitaion because pions are not annihilated finally. Just possibility of 
this process is recorded. 

Results obtained by the above described two independent methods of dilep- 
ton production simulation were nearly identical in all features (invariant mass 
distribution, pt distribution, azimuthal and rapidity distribution) tested so far. 



3 Invariant Mass Distribution 



Invariant mass distribution of dilepton pairs produced in ultrarelativistic heavy 
ion collisions can be divided into separate regions. Dileptons with invariant 
mass M > 3 GeV are produced mainly in the process of parton-parton annihi- 
lation. Region of M below 300 MeV is dominated by Dalitz decays of 7t^,t],uj 
and 7]' mesons. Also electromagnetic bremsstrahlung contributes to this low 
mass region of dilelectron pairs. Dilepton pairs at region of M close to IGeV 
are decay products of p,uj and mesons. In the region between 280 MeV 
and 1 GeV thermal production from tt'^tt" annihilation contributes. This list 
of dilepton production channels would be incomplete if we did not mention 
dileptons from decays of J/ip, ip' ^ind couples of DD mesons. 

TT+vr^ annihilation channel of dilepton production plays important role in 
the case of ultrarelativistic heavy ion collisions due to high multiplicities of 
secondary produced pions (above 2000 in central Pb-|-Pb 160 GeV/n collisions). 
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Fig. 2 Invariant mass distributions of dileptons obtained by the simulation. 
In this section invariant mass spectra of dileptons originating from tt+tt" 
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annihilation generated by the rescattering simulation are presented. In the pa- 
per mainly influence of formation time parameter on the shape of dilepton 
spectrum was investigated. Here dependence of the dilepton mass spectrum 
on the initial state of the pion gas is studied. 

Initial momenta and positions of pions in events produced by cascade gen- 
erator [0 are uncorrelated. It means that momentum of the pion does not 
depend on coordinates of its creation e.g. in transversal plane. Invariant mass 
spectrum obtained for these x — p uncorrelated events is shown on Fig.2a). 

Thick line on Fig. 2a) is obtained for same values of Tf and tj parameters 
of the rescattering model but for x — p correlated events. These events were 
obtained by the following modification of events produced by cascade generator 



Tl[| : For each pion scalar product of its transversal momentum and position 
in transversal plane pt ■ Xt was evaluated. In the case of negative numerical 
value of the scalar product, position of the pion xt was replaced by the oposit 
position in transversal plane —Xt. In this way momentum-space correlation 
was generated in the initial state of pion gas. Momenta of pions (not changed) 
are oriented predominantely out from the center of transversal plane. 

Resulting invariant mass spectrum of dileptons calculated for x — p cor- 
related events differs significantly from the spectrum of uncorrelated events. 
Low mass dilepton pairs are strongly enhanced in the case of x — p correlated 
events. In Fig. 2b influence of formation time parameter r/ on invariant mass 
spectrum of dileptons is shown for x — p correlated events. Small values of 
formation time parameter enhance production of dilepton pairs in p-meson 
region as it was found and discussed paper [|1^ . 



What is physical origin of the low mass dilepton enhancement in the case 
of X — p correlated events ? Why non-zero values of formation time suppress 
dilepton pairs in p-mass region compared to the low mass region pairs ? 

Two independent effects originating from the finite size of pion gas can 
serve as explanation. First one was already discussed in the work |jl2|. This 
consideration is reproduced very shortly here. Formation time [|1^ excludes 
pions from the rescattering process and therefore also from the possibility of 
annihilation into dilepton pairs. Faster pions are forbidden to interact for 
longer time due to the Lorentz dilation of the formation time implemented in 
the rescattering program. Due to the finite size of the pion gas faster pions 
can escape from the interaction region without any rescattering or annihilation. 
Thus invariant mass spectrum if dileptons is enhanced at low mass region by 
interactions of slow pions. 

At the same time another mechanism also resulting from finite size of the 
expanding pion gas influences invariant mass spectrum of dileptons. Invariant 
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mass quantity depends not only on the total size of momenta of pions being 
scattered but also on the angle between the momenta: 



M = ^{E^ + E2Y-{p^+p^y = pE^E2 + 2ml-2p^p2Cos{(P) (4) 

For small values of scattering angle 0, it means for pions moving in parallel 
directions, invariant mass quantitiy (^) is small. At the same time for pions 
moving in oposit directions invariant mass is large for the same values of mo- 
menta of pions. This effect can play significant role if the mean free path of 
interacting particles is comparable with the initial size of the expanding gas. 
This seems to be the case of the pion gas created in ultrarelativistic heavy ion 
collisions. 

Momentum-space correlation leads to the higher probability of collisions 
with small scattering angles and therefore it enhances production of low mass 
dileptons (see Fig. 2a). In the case of ultrarelativistic heavy ion collisions 
momentum-space correlation at the initial stage of pion gas can be generated 
also by transversal flow in prehadronic phase. 

4 Rapidity Distribution 

Rapidity distribution of dilepton pairs produced by pion annihilation is deter- 
mined by rapidity distribution of momentum sum Ptttt = f>i +P2 of anninhilating 
pions {pi,P2 are momenta of tt"*" and 7r~). This quantity is conserved in the 
process of tt+tt" annihilation, p^^r = pee- In Fig. 3a) and b) rapidity distribu- 
tions of pions before and after the rescattering process are shown. Rapidity 
distribution of pions is not influenced significantly by the rescattering process. 
However resulting rapidity distribution of dilepton pairs exhibits unexpected 
minimum at midrapidity surrounded by two maxima. 



Pb+Pb 160GeV/n b=7tm Pb+Pb 1 60 GeV/n b=7fm Pb+Pb 1 60 GeV/n b=7fm 




0.5 Y/Yb„, 1.0 0.5 Y/Yb,,„ 1.0 0.5 Y/Yb„„ 1.0 



Fig.3 Rapidity distribution of dileptons c) and pions before/after rescattering. 
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The result shown on Fig.3 is obtained for the simulation of Pb+Pb 160 
GeV/n events with values of rescattering model parameters Tf = 0.5 fm/c, 
ti = 1.0 fm/c. In this case number of collisions per pion is close to 1.0 and 
minimum in the rapidity distribution of dileptons is strongest. For higher 
collision rates (obtained for smaller vaues of r/ and ti parameters) the minimum 
is significantly smaller. For number of collisions per pion above 10 coll/vr the 
minimum disappeares at all. Such behaviour indicates that the minimum in 
the rapidity distribution of dileptons is generated by non-equilibrium process 
in the expanding pion gas. This effect is studied in PhD work 



From experimental point of view data on rapidity distribution of dileptons 
produced in p — A or A — A collisions are rather rare. It is not excluded that 
the only rapidity distribution of dileptons which seem to originate from vr+vr" 
annihilation was obtained in pioneering experiments of DLS collaboration at 
Bevalac accelerator in Berkeley [p!5|. 



5 Azimuthal Asymmetry of Dileptons 

Azimuthal asymmetries in transverse momentum distributions of particles have 



been clearly identified in ultrarelativistic |16[ non-central heavy ion collisions. 
This phenomenon is well understood as a consequence of collective behaviour 
of nuclear matter or explained by absorption of secondary produced particles in 
the spectator matter Recently also azimuthal asymmetries in transverse 
momentum distributions of less abundant hadrons - K mesons and A baryons 



T8[] have been identified in HIC experiments. However azimuthal asymmetries 
in transversal momentum distribution of dileptons have been addressed neither 
experimentally nor theoretically so far. 
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Fig. 4 Asymmetry in transversal momentum distribution of dileptons. 
From theoretical point of view main mechanisms generating azimuthal 
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asymmetries in transversal momentum distributions of hadrons are not appli- 
cable for dileptons. Dileptons leave freely collision volume after being produced 
without any collective final state interactions or absorption process. 

In spite of this the computer simulation described shortly in Section 2 pre- 
dicts significant second-order asymmetry in transversal momentum distribu- 
tion of dilepton pairs. Result shown in Fig.4 was obtained from the simulation 
of 40 non-central Pb-|-Pb 160 GeV/n 6 = 7 fm events. For each produced dilep- 
ton pair orientation of momentum pee = Pe+ + Pe- in respect to the reaction 
plane (impact parameter) was determined. Histogram of azimuthal angles 0ee 
of dilepton pairs was filled and significant second order asymmetry was found. 
The fit of the histogram shown in Fig.4 to the function: 

= S'o[l + S2-cos(20)] (5) 

gave numerical value of the asymmetry coefficient 5*2 = 0.093 ± 0.004. Corre- 
sponding value of Rp parameter (simplified for our purposes) used to charac- 
terize squeeze-out effect of hadrons at Bevalac energies is: 



«,^§^ 1.202 (6) 

Orientation of the reaction plane is parallel to y axis in our simulation, 
azimuthal asymmetry of dileptons is oriented in the reaction plane. From the 
analysis of second order asymmetry in the transversal momentum distribution 



of pions it comes out ^ that the asymmetry is strongest for the rate of 
collisions per pion close to 4.0 coll/vr. For higher collision rates strength of 
the asymmetry decreases though it does not tend to vanish in thermal equi- 
librium limit. The same behaviour is expected for the azimuthal asymmetry 
of dileptons presented here for the first time. 

Mechanism of generation of the azimuthal asymmetry of dileptons is well 



understood at present. The asymmetry is closely related to the spatial 
distribution of the pion gas in transversal plane which is asymmetrical in non- 
central events. 



6 Summary 

Two independent methods of the simulation of dilepton production from the 
expanding pion gas were found to generate nearly identical results. Invariant 
mass spectrum of dileptons produced via vr+vr" annihilation channel was shown 
to be sensitive to momentum-space correlation in initial stage of the pion gas. 
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Generated rapidity distribution of dileptons exhibits unexpected minimum 
at midrapidity. The minimum seems to be a consequence of a non-equihbrium 
process in the expanding pion gas since it becomes weaker and disappears for 
higher number of colhsions (above 10 coU/vr) in the pion gas. 

For non-central colhsions second order asymmetry in azimuthal distribu- 
tion of dileptons is found. This asymmetry is a consequence of asymmetri- 
cal overlapping region of the colliding nuclei in transversal plane. Azimuthal 
asymmetry of dileptons might be accessible experimentally if the reaction plane 
orientation is determined from hadronic signal on dilepton detector. 

7 Conclusions 

In this contribution we study expansion of interacting pion gas created in ultra- 
relativistic heavy ion collisions. Non-equilibrium and asymmetrical features of 
the expansion exhibit itself in rapidity and azimuthal distributions of particles 
participating or created in the rescattering process. Rescattering phenomenon 
is not related only to the hadronic phase of nuclear matter produced in HIC. 



It can equally be present also in prehadronic phase where it can lead 
to similar effects as those revealed in the pion gas simulation. Investigation 
of the non-equilibrium expansion of the parton gas possibly created also in 
non-central collisions might shed light on other interesting aspects |^ of HIC 
research. 



Content of this article was presented at HEP'97 conference in Jerusalem [E3 
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